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Abstract. The influence of both major and minor mergers is expected to significantly affect 
gradients of stellar ages and metallicities in the outskirts of galaxies. Measurements of observed 
gradients are beginning to reach large radii in galaxies, but a theoretical framework for con¬ 
necting the findings to a picture of galactic build-up is still in its infancy. We analyze stellar 
populations of a statistically representative sample of quiescent galaxies over a wide mass range 
from the Illustris simulation. We measure metallicity and age profiles in the stellar halos of quies¬ 
cent Illustris galaxies ranging in stellar mass from 10^° to accounting for observational 

projection and luminosity-weighting effects. We And wide variance in stellar population gradi¬ 
ents between galaxies of similar mass, with typical gradients agreeing with observed galaxies. 
We show that, at fixed mass, the fraction of stars born in-situ within galaxies is correlated with 
the metallicity gradient in the halo, confirming that stellar halos contain unique information 
about the build-up and merger histories of galaxies. 


1. Introduction 

Stellar halos are diffuse regions of stars ubiquitously found surrounding galaxies (with 


a few notable exceptions. 

van Dokkum et al. 

2014 

). They are observed to extend to 

many times a galaxy’s effective radius ( 

Martmez-Delgado et al. 

2010 

), where dynami- 


cal timescales are very long compared to the ages of their host galaxies (Eggen et al. 


19621. Detections of streams and tidal features in integrated light images of the halos in 
the Milky Way (Helmi et al.|1999), Andromeda (Ibata et al.|2001 ), and external galaxies 


as galaxies grow hierarchically, a picture anticipated by numerical simulations (e.g., 

John- 

ston et al.p996 

Bullock &: Johnston|2005 

Cooper et al.|2010 

Pillepich et al.|2014 

). The 


long timescales at large radius can help preserve the information content of each galaxy’s 
merger history, which has been used as motivation for observations of kinematics and 
populations of the stellar halo (Bell et ^|2008 Schlaufman et al(||2009 1. 


Stellar populations - characterized by the metallicities and ages of stars - are clues 


to when and in what systems stars originally formed (Greene et al. 2015). The mass- 


metallicity relation in galaxies (Tremonti et al. 2004), combined with hierarchical ac¬ 
cretion, implies that metallicity gradients should indicate the relative contributions of 
more-or-less massive progenitor systems to particular regions of a galaxy and its halo. 
Different merger histories should leave distinct imprints on observed stellar population 
gradients ( Hirschmann et al.|2015 ), although stars formed in-situ, or within the galaxy in 
which they now reside, may still be significant in the inner regions of the halo ( [Font et al. 
2011). For massive, early-type galaxies (usually red, quiescent, and elliptical), the two- 


phase formation scenario outlines the contributions of mergers and in-situ star formation 


to the formation of a metallicity gradient (Spolaor et al. 2010 Pastorello et al. 2014). 


The early phase of dissipative collapse leads to an steep, negative gradient from in-situ 
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stars, while the later phase involving accretion of stars from smaller satellite galaxies 
tends to flatten these gradients ( Kobayashi]|2004 1. 

Observations of stellar population gradients are beginning to accumulate for large sam¬ 
ples of galaxies ( Spolaor et al.|[20T0 Pastorello et al.||2014 1 and are reaching continually 
larger radii through integral-field spectroscopy ( Delgado et al.|20l5 Greene et al.||2015 ). 
Observations find a wide variety of gradients between galaxies with similar masses and 
morphologies, which is to be expected if particular merger histories shape the gradients 
in a stochastic way. Unfortunately, the observational need for detailed predictions con¬ 
necting gradients to merger histories have so far outpaced numerical simulations. Due to 
limitations in computing resources, hydrodynamical simulations have so far only resolved 
stellar halo populations around individual galaxies ( Abadi et al.|2006 Hirschmann et al. 
2015 Cooper et al.|2015 ), a method which lacks the statistical power required to replicate 
the observed halo-to-halo variations. Large samples of galaxy halos have been produced 
in N-body simulations, some of which include semi-analytic models and stellar-tagging 
techniques ( Cooper et al.|2010| , but these do not include baryonic physics effects, which 
may shape the dark matter distribution, and must rely on complicated fitting functions 
to generate “stellar particles” with realistic orbital properties ( Bailin et al.||2014 ). 

In this work, we measure stellar population gradients in a sample of quiescent galaxies 
from the hydrodynamical cosmological simulation Illustris. With the statistical power 
of Illustris’ large cosmological volume and its self-consistent model for star and galaxy 
formation, we are able to match the observed galaxy-to-galaxy variance in gradients and 
show that these variations are indeed reflective of different galactic merger histories. 


2. Methods 


2.1. The Simulations 


The Illustris simulations (Vogelsberger et al. 2014a 


2014a 

Genel et al. 

2014| 

Nelson et al. 


20151 are a suite of N-body-l-hydrodynamical cosmological simulations (106.5 Mpc on a 


side), run at multiple resolutions with the adaptive mesh code AREPO (Springel 2010 
Vogelsberger et al.|2013 1. The simulations model key physical processes for the formation 
of galaxies, including stellar formation, evolution, and feedback, chemical enrichment, 
radiative cooling, supermassive black hole growth, and feedback from AGN. The highest 
resolution run (Illustris-1, hereafter simply Illustris) has a mass resolution of moM = 
6.26 X 10® Mq and mharyon 1-26 x 10® Mq for the dark matter and baryonic components, 
respectively. At 2 = 0 gravitational forces for stellar particles are resolved to a softening 
length of 0.7 kpc. Illustris was run from z = 127 to z = 0 using ACDM cosmological 
paramters consistent with WMAP9 (Da = 0.7274, = 0.2726, h = 0.704, Hinshaw 

eFalllM^ . 


At z = 0, the Illustris volume contains more than 4 x 10"^ well-resolved galaxies (Vo- 
gelsberger et al.|2014a l, with a reasonable diversity of morphologies and colors, including 
early-type and late-type galaxies ( Torrey et al.|2015 |. The most massive central galaxies 
(as identified by the FOP and SUBFIND algorithms, Springel et al. 2001 Dolag et al. 
20091 reproduced in the simulation have stellar masses within their stellar half-mass radii 
of M* ~ 1 — 2 X IO^^Mq. The simulation reproduces the observed z = 0 mass-metallicity 
relation in galaxies ( Vogelsberger et al.|2014b ), and a reasonable relation between galaxy 
mass and stellar ages once luminosity weighting is properly taken into account. 

2.2. The Quiescent Galaxy Sample 

Our goal to reproduce measurements in the outer regions of early-type galaxies moti¬ 
vates the selection criteria of our sample. We select central galaxies (which are not satel- 
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lites/subhalos of more massive parents) with stellar masses M* ^ Each galaxy 

in this mass range is resolved with at least a few hundred star particles beyond 8 effective 
radii, ensuring that we constrain gradients all the way throughout the low-density outer 
regions and up to 10 times the effective radius. 

In this work, we examine the properties of simulated quiescent galaxies, in comparison 
to observed early-type galaxies. Of all Illustris galaxies with M* > 10 ^°Mq, we select 
the 352 quiescent galaxies with specific star formation rates SSFR ^ 10“^^'® yr“^ within 
twice the stellar half-mass radius. See Cook et al. (In Prep.) for details. 

2.3. Fitting Stellar Population Gradients 

We simulate observational projection effects by projecting the radii of star particles from 
the center of their host galaxies against a random line-of-sight. We also account for 
observational biases by weighting the contribution of each star particle relative to its 
V-band lumionsitjEl Star particles identified by SUBFIND as bound to smaller subhalos 
or satellites are not included in the analysis. 

Using these observational considerations, we measure the 2-D azimuthally-averaged 
values of stellar metallicity and age in five logarithmically-spaced radial bins over a 
chosen projected radius range. We focus on three particular ranges in units of the V- 
band effective radius (i?el31), which we label: the inner galaxy (O.I - 1 Re), the outer 
galaxy (1-2 Re), and the stellar halo (2 - 10 Re). We then calculate a logarithmic 
gradient by fitting a line to these averages, equally weighting each radius bin: 

f{r) = f{Re)+yf\og,o{r/Re), ( 2 . 1 ) 

with / the mean log-metallicity ([Z/H]) or the mean age (in Gyrs) in each bin. 


3. Results and Implications 

3.1. Comparisons to Observed Gradients 

The measured metallicity and age gradients in the inner and outer galaxy ranges (< 
2Re) are shown as a function of central velocity dispersion (erg, within |i?e) in Fig. 1 
Observations can constrain gradients in these radius ranges for large samples of galaxies 


( Spolaor et al.|2010 Pastorello et al.|20T4 Greene et al. 2015 Delgado et al.|201^ , but so 

far there are only a few individual cases for comparison in the stellar halo (> 2Re) range. 
We compare our measurements to observations of individual galaxies from [Spolaor et al.| 
(2010) (inner galaxy) and stacked measurements by Greene et al. (2015) (outer galaxy). 


The typical values of our metallicity and age gradients agree with observations: metal¬ 
licity gradients are negative - outer regions are more metal poor - while age gradients 
are roughly flat - inner and outer regions have similar average ages. Of particular impor¬ 
tance is the fact that our measurements also reproduce the observed scatter in gradients 
between galaxies of similar masses. Illustris, with its large statistical sampling of galax¬ 
ies, is uniquely able to replicate this galaxy-galaxy variance, whereas more individualized 
simulations cannot replicate the wide variety of observed gradients. 

f Throughout the paper, M* denotes the total mass in stars within the stellar half-mass 
radius (Rhaif) after removing the contributions of gravitationally-bound satellites identified by 
SUB FIND. 

j: Star particles in each galaxy are assigned luminosi ties in several com mon observational 
bands using single-age stellar population SED templates (Torrey et al.|2015l. 

^ Each galaxy’s effective radius is the projected radius within which one-half of its total 
V-band luminosity is located. This is calculated as an average over 100 random lines of sight. 
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Figure 1. Measurements of metallicity (top) and age {bottom) gradients in Illustris quiescent 
galaxies (black cross marks). These are calculated in the inner (0.1 - 1 Re, left) and outer (1 - 
2 Re, right) galaxy ranges, and include projectio n along a random lin e of sight and luminosi- 
ty-weighting. We compare to the measurements of Spolaor et al. (20101 and Greene et al. ( |2015 1 
in the inner and outer galaxy, respectively. 


3.2. Relating Metallicity Gradients to Merger Properties 

Here, we study whether the scatter in stellar population gradients among Illustris galaxies 
at fixed mass - which matches the scatter observed in early-type galaxies - can be 
connected to the particular merger histories of each galaxy. We quantify the contribution 
of mergers to a galaxy’s build-up in terms of the in-situ fraction: the mass fraction of 
a galaxy’s stars which were formed within the galaxy or its main progenitor branch. 
Galaxies with low in-situ fractions were primarily built-up from mergers and accretion of 
smaller systems, while galaxies with high in-situ fractions have had little influence from 
mergers. 

In the hierarchical model of galaxy formation, galaxies grow continually through merg¬ 
ers with increasingly massive neighbors. This effect, reproduced in many cosmological 
simulations and seen clearly in Illustris ( Rodriguez-Gomez et al.|2015 ), results in a strong 
trend towards lower in-situ fractions in more massive galaxies. To disentangle this effect 
from the individual histories, we can focus our analysis to small ranges in stellar mass. 
Over a sufficiently narrow window, mass differences are small enough that the in-situ 
fraction should be driven primarily by the particular merger history of each galaxy. 

In Fig. we show the metallicity gradient in the stellar halo (2 - 10 Re). Galaxies 
have overall negative metallicity gradients and there is a large scatter at fixed mass, just 
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2 - 10 Re 




Figure 2. Top: The metallicity gradient in the stellar halo (2 - 10 Re) of the quiescent Illustris 
sample, versus galaxy stellar mass. As at smaller radii, the metallicity gradient shows wide 
galaxy-galaxy variance within narrow mass ranges, such as between 10^° ® and (shaded 

region). Bottom: The correlation between galaxy in-situ fraction and stellar halo metallicity 
gradient, for galaxies between in the mass range given above. Galaxies with lower in-situ fractions 
(mergers more significant in growth) have flatter metallicity gradients. 


as in the inner and outer galaxy regions. When selecting galaxies in the narrow mass 
window from 10^° ® to 10 ^^Mq and comparing to the in-situ fraction, we see that there 
is a correlation. Low in-situ fractions are associated with flatter metallicity gradients in 
the stellar halo, while galaxies with higher in-situ fractions have steeper halo gradients. 

This result suggests that much of the scatter in observed metallicity gradients can be 
explained through a variety of galactic merger histories. The process of building-up a 
galaxy through mergers appears to flatten the metallicity gradient. Measuring gradients 
into the outer stellar halo should provide clues about the relative importance of mergers 
in a particular early-type galaxy’s evolution (see also, Pillepich et al.|[20l4 ). 

We note that this correlation is significant only in the stellar halo regions (2 - 10 Re) 
of our simulated galaxies. Metallicity gradients in the inner and outer galaxy regions 
(< 2Re) show no association with the in-situ fraction. This emphasizes the importance 
of measuring gradients into the outer stellar halo (well beyond 2Re), where the signifi¬ 
cantly longer dynamical timescales result in more lasting imprints from mergers. Stellar 
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population gradients in the halo retain the information content of merger histories; it is 
still unclear what leads to the wide scatter in gradients in the interior regions. 
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